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On the basis of the models of various developed high-temperature fuel-cell heat-engine 
hybrid systems, a unified model of hybrid systems is proposed. General expressions for 
the power output and efficiency of hybrid systems, high-temperature fuel cells such as 
solid oxide fuel cells (SOFCs) and molten carbonate fuel cells (MCFCs), and heat engines 
including the Brayton, Otto, Diesel, Atkinson, Braysson, and Carnot engines are, respec¬ 
tively, derived by using the theories of electrochemistry and non-equilibrium thermody¬ 
namics. The effects of main irreversible losses existing in real fuel cells and heat engines 
on the performance of hybrid systems are investigated. The general performance char¬ 
acteristics and optimal operating regions of some of the key parameters of hybrid systems 
are discussed in detail. A variety of special typical cases are discussed. The important re¬ 
sults in the literature can be readily reproduced, and the interesting findings of our study 
are presented. 

Copyright © 2013, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights 

reserved. 


1. Introduction 

In general, fuel cells can yield the highest efficiency and 
lowest emission for any known fossil-fueled power generation 
technology [1], Consequently, they have received considerable 
attention in the world. Especially for high-temperature fuel 
cells, such as solid oxide fuel cells (SOFCs) and molten car¬ 
bonate fuel cells (MCFCs) [2,3], they boast advantages of coal 
and biomass gasification [1,4,5], fuel internal reforming [6,7], 
fuel flexibility [8,9], C0 2 capturing [10;, high-quality heat 
coproduction [11], etc. However, high-temperature fuel cells 


face some market entry challenges, ranging from the cost to 
the durability [12-14], Researchers are making an effort in the 
development of suitable materials and fabrication of SOFCs 
[15-17], so that SOFCs could be suitable for small-scale resi¬ 
dential market applications at sufficiently low cost of $1000/ 
kW [18]. A 3.5-kW SOFC cogeneration system for residential 
and commercial applications with an electrical efficiency of 
about 50% and overall efficiency of approximately 80% has 
been developed by Hydrovolt [19]. The cost of MCFCs is 
generally low in comparison with that of SOFCs when cheaper 
materials and inexpensive fabrication techniques are used. 
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The power density of MCFC systems operating at 650 °C, 
however, is relatively low, and is approximately 1.5 kW m~ 2 
[20]. It follows that MCFC will likely be commercialized at sizes 
greater than 100 kW, and such commercialization will 
certainly increase the cost. Hence a more appropriate strategy 
is to develop smaller MCFC systems for distributed combined 
heat and power (CHP) applications with a typical electrical 
power output between 250 kW and 1 MW [20], The cost of fuel 
cell systems can be reduced not only by better-quality mate¬ 
rials [21—23] but also by further-advanced technologies [24] 
that include the integration of high-temperature fuel cells 
with heat engines. 

For purposes of low pollutant emission and high fuel 
economy, a hybrid system including fuel cells and Diesel or 
Otto engines was proposed in 1999 [25 . Since up to a half of 
the fuel energy of high-temperature fuel cells is waste heat, 
numerous investigators have carried out the research on high- 
temperature fuel-cell based hybrid systems in order to in¬ 
crease the overall efficiencies of fuel cells. Research activities 
are exemplified by (a) using high-temperature fuel cells as 
topping units within gas turbine cycles [26,27], (b) overcoming 
lower performances of high-temperature fuel cells than those 
of low-temperature ones by constructing fuel cell-Carnot heat 
engine hybrid systems '28], (c) recovering waste heat from the 
cell by the design of SOFC-Stirling or MCFC-Stirling hybrid 
systems [29], and (d) proposing a Brayton-cycle thermal- 
management system to reduce the hydrogen fuel flow for 
cooling [30], High-temperature fuel-cell based hybrid systems 
also include CHP systems, combined cooling, heat, and power 
(CCHP) systems [31,32], SOFC-proton exchanger membrane 
fuel cell (PEMFC) systems [33], and fuel-cell and Rankine or 
Kalina cycle hybrid systems [34-36]. 

The integration of heat engines and fuel cells are capable of 
not only recovering energy from fuel-cell exhaust gases, but 
also generating extra power. The efficiencies of SOFC- or MCFC- 
based hybrid systems can attain up to 70%-90% [37]. By 
comparing the technologies of fuel cells and internal com¬ 
bustion engines in vehicles, the former enjoys simpler designs, 
higher reliabilities, noiseless operations, higher efficiencies, 
less environmental impact, and lower energy consumption 
[38-40], However, fuel cells may not supply sufficient energy 
during peak power demand periods and transient events, and 
cannot save power and recycle braking energy, either. The 
optimal combination of the fuel cell, internal combustion en¬ 
gine, and electric energy storage and conversion devices is 
proposed. It is more flexible with its efficiency being less load- 
dependent than individual units un-combined [25,41,42]. 

In the present paper, the general performance character¬ 
istics of high-temperature fuel cells, heat engines, and their 
hybrid systems are described. The advantages of hybrid sys¬ 
tems are demonstrated in numerical results, with contents 
organized as follows. In Section 2, a unified model of high- 
temperature fuel-cell heat-engine hybrid systems is estab¬ 
lished. The power output and efficiency of hybrid systems are 
derived. In Section 3, some general performance characteris¬ 
tics of hybrid systems are revealed and their parametric 
optimal criteria are determined. Many interesting cases are 
discussed in detail in Section 4. The performance character¬ 
istics of several developed hybrid systems are presented. 
Finally, some important conclusions are drawn. 


2. A unified model of high-temperature fuel¬ 
cell heat-engine hybrid systems 

A generic hybrid system primarily consists of a high- 
temperature fuel cell and a heat engine, as shown in Fig. 1 
[28,29], where the high-temperature waste heat generated in 
the fuel cell is utilized as the heat input of the heat engine, the 
waste gases produced in the fuel cell are used to preheat the 
fuel and the oxidant through a heat exchanger, P e and P h are 
the power outputs of the fuel cell and the heat engine, 
respectively, q-i and q 2 are the rates of heat flow from the fuel 
cell at temperature T to the working substance of the heat 
engine and from the working substance of the heat engine to 
the environment at temperature T 0 , respectively, and qi is the 
heat leak rate from the fuel cell to the environment. Below, we 
first introduce the developed models of high-temperature fuel 
cells and heat engines, and then establish a unified model of 
high-temperature fuel-cell heat-engine hybrid systems. 

2.1. Power output and efficiency of high-temperature 
fuel cells 

When studying the performance of high-temperature fuel 
cells, one has to consider the thermodynamic and electro¬ 
chemical irreversibilities which result from the activation 
overpotential (rj act ), ohm overpotential (rjohm), and concentra¬ 
tion overpotential (rj con ) [43,44]. A general equation 

Pact T Vohm T Peon _f j ^ 0 T jA 

RT In / / PH 2 PH 2 Q.bPro 2 .bPcQ 2 

"eF VPH 2 ,bPH 2 0Pco 2 PcO 2 .b V Po 2 .b J 

(1) 

can summarize the three overpotentials of high-temperature 
fuel cells based on the theoretical and experimental results, 
where i and i 0 are, respectively, the current density and the 
exchange current density, pj and Sj are the specific resistivity 
and the flowing length in layer j, respectively, R is the uni¬ 
versal gas constant, n e is the number of electrons transferred 
in reaction, F is Faraday’s constant, p;(i = H 2 , H 2 0, 0 2 ,...) is the 
partial pressure of component i, and p j|b is the partial pressure 
of component i in bulk phase on the surface of electrodes. 





Fig. 1 — A schematic diagram of the unified model of high- 
temperature fuel-cell heat-enigne hybrid systems. 
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Equation (1) can be expressed in several forms for differently 
working conditions, which can be seen from the following 
discussion for three overpotentials. 


2.1.1. Activation ouerpotential 

Activation overpotential occurs due to the energy barrier 
which must be overcome by reactants to allow the electro¬ 
chemical reaction to happen. An accurate relation [45—48] 


i 

— = exp 
lo 


(an e F \ 


V RT Vact ) 

- exp 


(a - l)n e F 
RT 


Pact 


( 2 ) 


between the activation overpotential and the current density 
is derived from Butler-Volmer equation [49]. At high- 
activation polarization, Eq. (2) will be simplified into the 
Tafel expression [50,51] 
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RT 
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In 



(3) 


Pi 


— ajexp 



(9) 


while there are other expressions for the electrolyte’ specific 
resistivity [59,60] 

p e = T e a e expf!^) (10) 

or [61] 


p e = a e exp 





( 11 ) 


The specific resistivity of interconnectors can also be a 
constant to fit experiment data [29]. In Eqs. (9)-(ll), aj, b,, a e , 
and b e are constants. If there is not a generally applicable 
mechanistic equation to calculate the total internal resistance 
caused by ohmic losses, the ohm overpotential can also be 
expressed as [62] 


while at low-activation polarization, Eq. (2) will then be 
simplified into the well-known linear potential-current rela¬ 
tion [52,53] 


RT i 


(4) 


Another linear relationship between the activation over¬ 
potential and the current density obtained by experimental 
approach is [54] 


Vact = aiexp 



(5) 


In Eqs. (2)— (5), a is the charge transfer coefficient of elec¬ 
trodes, a and b are two constants, and i 0 can be derived by a 
semi-empirical method and is given by [45-48,52,53] 


• ( \ ai ( \ ° 2 ( \ Cl / \ fl3+C 2 

1o = 7o(Ph 2 (co>J (Ph 2 o) [Po 2 ) (.Pco 2 J exp 


Eact\ 

RT J 


( 6 ) 


or [55] 


'Eact /A _iy 
R \To Tjy 

(7) 

where ci = c 2 = 0 for the anode, a t = a 2 = a 3 =0 for the cathode, 
7 o is the coefficient of electrodes, E act is the activation energy 
of electrodes, T 0 is the reference temperature of high- 
temperature fuel cells, and p H2 (co), Po 2 , Ph 2 o, and p C o 2 are the 
partial pressures of H 2 (CO), 0 2 , H 2 0, and C0 2 , respectively. 

2.1.2. Ohm ouerpotential 

Ohm overpotential occurs due to the resistance to the flow of 
ions in the electrolyte and the resistance to the flow of elec¬ 
trons through electrodes and interconnectors. By Ohm’s law, 
the ohm overpotential existing in the anode, cathode, elec¬ 
trolyte, and interconnectors of high-temperature fuel cells is 
[47,56] 


• ( \ Ql / \ a 2 / \ C 1 / \ a 3+ c 2 

io = To(Ph 2 ( co)J [Ph 2 oJ [Po 2 J [Pco 2 J exp 


Vohm = iJ2 p i 6 i’ ( 8 ) 

i 

where pj is dependent on the temperature Tj of layer j, which is 

[47,57,58] 


^ + t 2 T + y 3 i + y 4 Ti + 7 5 T 2 + y 6 i 2 , (12) 

where y lt y 2 , ..., 76 are constants determined by experiment 
data. 


2.1.3. Concentration ouerpotential 

Concentration overpotential occurs when the diffusion rates 
of reactants and products through electrodes are slower than 
the electrochemical reaction [51]. For MCFCs, the reaction 
takes place at low polarization so that the concentration loss 
is negligible, compared with other losses 44,52], The diffusion 
transport in a porous material can be described by Fick’s 
model, dusty gas model, or Stefan—Maxwell model [63,64]. 
Fick’s model, which is simple and is used more frequently, is 
adopted here. According to Fick’s model, the concentration 
overpotential can be expressed as [46,65] 


Vcon " 


n e F [ ^ n e FD H2 p H2 ,b J / \ 


RT 
" 2n„F 


■In 


Pc a 
Po 2l b 


Pea 

Po 2 ,b 


1 exp 


RTd c i 

2n c FDo 2 p C a 


(13) 


where 5 a and <5 C are, respectively, the thickness of anode and 
cathode electrodes, p ca is the pressure of the cathode elec¬ 
trode, and Dh 2 , Do 2 , and D H2 o are, respectively, the effective 
diffusion coefficients of H 2 , 0 2 , and H 2 0, which depend on the 
temperature and the structure of fuel cells. 

In some literature, Eq. (13) was simplified into [51] 


Vcon 


RT, 
—=ln 
n„F 



RT 

2t%F‘ 


In 1 



or [42,45,56,57] 


(14) 




(15) 


where i L-H2 and i L02 are, respectively, the limiting current 
densities at which the hydrogen and the oxygen are used up 
at a rate equal to their maximum supply speeds, and i L is the 
limiting current density at these electrodes. There is another 
common expression to model the concentration over¬ 
potential [66], i.e. 
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Vcon = mexp(ni), 


(16) 


where m and n are constants. 

The open-circuit voltage of high-temperature fuel cells is 

[51,67,68] 


—Ag°(T) R T ln ( PHiPofPco^ 

n e F n e F ^ Ph 2 oPco 2 J 


(17) 


where Ag°(T) is the molar Gibbs free energy change of the 
chemical reaction at temperature T and standard pressure, 
and p c c a 0 and pff 0 are the partial pressures of C0 2 at different 
electrodes, respectively. Considering the influence of irre¬ 
versible losses above on the performance of high-temperature 
fuel cells, one can derive the power output and efficiency of 
high-temperature fuel cells as 


P e = iA c (E - V aa - V ohm - V con ) 


= iA c 


-Ag°(T) RT t PH 2 Po / >g 1 o ; \ 
n e F n e F ^ Ph 2 oPco 2 J 



•V- s RT 1 


n P F 


PH 2 PH 2 o,bPro 2 .bPco 2 

PH 2 .bPH 2 oPS 2 Pco 2 ,b 



(18) 


and 


Pe 

Ve -iA c Ah°(T)/(n e F) 

n e F I" —Ag°(T) RT / WffPcoA fi\ 

—Ah°(T) [ n e F n e F ^ Ph 2 oPS 2 J J U J 

-y 1 RP j n ( Pu 2 PH 2 o,bPco 2 ,bPco 2 jPo 2 A 

nj n (^PH 2 ,bPH 2 oPc0 2 Pc0 2 ,b V PojJ J ’ 

(19) 

respectively, where A c is the surface area of bipolar plates and 
Ah°(T) is the molar enthalpy change of the chemical reaction 
at temperature T and standard pressure. 


2.2. Two modes of heat engine cycles 

A survey of the developed heat engines has disclosed that 
there mainly exist two classes of heat engine cycles. The first 
class of cycles are composed of two isothermal and two pol¬ 
ytropic processes such as the Carnot cycle consisting of two 
isothermal and two adiabatic processes [69], the Stirling cycle 
consisting of two isothermal and two isochoric processes [70], 
the Ericsson cycle which consists of two isothermal and two 
isobaric processes [71], etc. The performance of the hybrid 
systems composed of such a class of heat engines and high- 
temperature fuel cell has been analyzed [26,28,29,72]. The 
second class of cycles are composed of two polytropic and two 
adiabatic processes such as the Carnot cycle consisting of two 
isothermal and two adiabatic processes, the Diesel cycle 
consisting of an isobaric, an isochoric and two adiabatic pro¬ 
cesses [73,74], the Atkinson cycle consisting of an isochoric, an 
isobaric and two adiabatic processes [75,76], the Otto cycle 
consisting of two isochoric and two adiabatic processes [77], 
the Brayton cycle consisting of two isobaric and two adiabatic 
processes [78], the Braysson cycle consisting of an isobaric, an 
isothermal, and two adiabatic processes [79], etc. The perfor¬ 
mance of some of the hybrid systems composed of such a 


class of heat engines and high-temperature fuel cells has been 
also analyzed [29,34,80]. 

On the basis of the cycle models of heat engines afore¬ 
mentioned, we can establish a generic cycle model of heat 
engines which include not only the second class of cycles 
mentioned above but also some new cycle modes, as show in 
Fig. 2, where qi and q 2 are the rates of heat flow absorbed from 
the heat reservoir at temperature T and released to the cold 
reservoir at temperature T 0 by the working substance during 
two polytropic processes 3-4 and 6-1, respectively, 1—2S and 
4-5S are two reversible adiabatic processes, 1—2 and 4-5 are 
two irreversible adiabatic processes, 2S—4 and 5-1 are two 
polytropic processes, q R is the rate of heat flow from regen¬ 
erative process 5-6 to 2—3, and T f (i = 1, 2, 3, 4, 5, 6) are the 
temperatures of the working substance at state points i. As¬ 
sume that heat capacities of the working substance in two 
polytropic processes are, respectively, c H and c L . If c H = c p and 
Cl = c„, the cycle is the Diesel cycle, where c p and c„ are the heat 
capacities of the working substance at constant pressure and 
volume, respectively; if c H = c L = c p , the cycle is the Brayton 
cycle; if c H = c L = c v , the cycle is the Otto cycle; if c H = c L -> , 

the cycle is the Carnot cycle; if c H = c p and c L -> <*>, the cycle is 
the Braysson cycle; if c H —°° and c L = c p , the cycle is a new 
cycle [81]; if c H = c„ and c L = c p , the cycle is the Atkinson cycle 
[82], In addition, such a cycle model also includes the cycle 
mode of the turbine engine [83], which is often adopted in the 
high-temperature fuel-cell turbine-engine hybrid systems 
[24,26,57,84]. When the cycle is endoreversible, the perfor¬ 
mance of the Carnot cycle is the same as that of the Stirling 
and Ericsson cycles with the ideal regeneration [72]. It implies 
that the results obtained from the generic cycle model may be 
directly used to discuss the performance of the endoreversible 
Stirling and Ericsson cycles with the ideal regeneration. 

When irreversibilities from non-ideal adiabatic and regen¬ 
erative processes in the cycle are considered, one may obtain [85] 

mc H ln < mc L ln (20) 

13 11 

according to the second law of thermodynamics. So a 
parameter I may be introduced to describe the internal 



Entropy 


Fig. 2 - Temperature-entropy diagram of a generic cycle 
model of the heat engines consisting of two poly tropic and 
two adiabatic processes. 
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irreversible losses of the working substance [86] and 
0 < I = (T 1 /T 6 )(T 4 /T 3 ) Ch/Ci < 1. When I = 1, the cycle is endor- 
eversible. Some authors [87,88] have also defined the effi¬ 
ciencies of the compressor, turbine, and regenerator, 
respectively, to describe the internal irreversible losses of the 
cycle. If the efficiencies of these three processes are all equal 
to 1, the cycle is endoreversible. 

Assume that the heat transfer between the working sub¬ 
stance and heat reservoirs obeys the Newtonian law [89,90], 
one can obtain the heat flow rates in the cycle as [73,81] 


K h A h (T 4 - T 3 ) 
qi ln[(T-T 3 )/(T-T 4 )] 


fiic H (T 4 — T 3 ), 


( 21 ) 


K l A l (T 6 - TO 

q2 ln[(T 6 - T 0 )/(Ti - To)] 


mc L (T 6 -Ti), 


and 


( 22 ) 


_ K R A R [(T 5 - T 3 ) - (T 6 - T 2 )] 
qR ln[(T 5 - T 3 )/(T 6 - T 2 )] 


rh c h (T 3 - T 2 ) = mc L (T 5 - T s ), 
(23) 


where K H , K L , and K R are heat transfer coefficients in processes 
3-4, 6—1, and 5—6, respectively, and A H , A Ll and A R are cor¬ 
responding heat transfer areas. Combining Eqs. (21)-(23), one 
can obtain the efficiency and rewrite the rate of the heat input 
of the heat engine as 


T 1 [(T 4 /T 3 ) T /I — 1] 

r( t 4 - t 3 ) 


and 


(24) 


2.3. Power output and efficiency of the unified model 

As shown in Fig. 1, one part of the waste heat from the high- 
temperature fuel cell is released into the environment 
directly through conductive and/or convective heat transfer, 
and the other part is transferred into the heat engine. The heat 
leak rate may be expressed as [91]: 

q, = K,A,(T-To), (28) 

where Ki and Aj are, respectively, the heat leakage coefficient 
and effective area of the fuel cell. Combining the first law of 
thermodynamics and Fig. 1, one can obtain the heat input rate 
of the heat engine as [92] 


qi = ( - iA c Ah°(T)) j (n e F) - P e - q,. 


(29) 


Combining Eqs. (19) and (28), one can rewrite Eq. (29) as 

qi/(AK H ) = bi[(l - Vefi ~ b 2 (T - T 0 )], (30) 

where b 4 = -A c Ah°(T)/(n e FAK H ) and b 2 = KiAin e F/(-A c Ah°(T)). 

Using Eq. (30), one can obtain the power output and effi¬ 
ciency of the hybrid system as 

} (31) 

and 


„ „ „ iA c Ah u (T) f 

P = P e + P h = - _ n p + Vh.r 


1 ~ Ve ~ 


b 2 (T - T 0 ) 


V = P/ 


( iA c Ah°(T)\ 


= Ve + Vh,m 1 


Ve 


b 2 (T - To) ~ 
i 


(32) 


qi = 


AK h (T 4 - T 3 ) 


T - T 3 K h t(T 3 - T t /x) + T 1 (T 4 /T 3 ) t /I - T, K h (T 1 (T 4 /T 3 ) t )/I - T 0 ’ 

T-T 4 K R ( 7 -1) (T 1 (T 4 /T 3 ) t )/I — T 4 /x 7 K l Ti-To 


(25) 


respectively, where A = A H + A L + A R is the total effective heat 
transfer area of the heat engine, y = c H /c L , and x = T 4 /T 2 . 

In order to derive the maximum efficiency for a given rate 
of the heat input, we introduce a Lagrangian function as 

L = rj h + 7qi, (26) 


respectively. For SOFC-based hybrid systems, )j e in Eqs. (31) 
and (32) may be expressed as [48,93] 


Ve 


-RT —Ag°(T) J ph 2 VP^ \ 
Ab»(T) RT ^ P h 2 o J 


2sin h 1 



where 7 is the Lagrangian multiplier. Using the extreme 
conditions 


' 9L \ 


= 
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Ti,t 3 
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+ In 


in e F 

Rf 


5 i a i ex P 




+ l l n 



(33) 


and for MCFC-based hybrid systems, rj e in Eqs. (31) and (32) 
may be expressed as [94—96] 


and Eq. (25), one can obtain the optimum values Ti^, T 3iV , and 
T 4j , of Ti, T 3 , and T 4 , respectively, for given q 1 /(AK H ), y, I, x, K H / 
K l , and K h /K r . Substituting T l 7J , T 3 7J , and T 4 ,, into Eq. (24), one 
can obtain the maximum efficiency r/ h m of the heat engine for 
a given qi/(AK H ). 


-RT [ —Ag°(T) f pH 2 ^P^Pco 2 ,ca ) 
Ve Ah°(T)\ RT + y P„ 2 oPc0 2 ,an J 


i i 

Io,a '0,c 


in e F 
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i5 e a e exp 
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( 34 ) 
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3. General performance characteristics and 
parametric optimum criteria 

Numerical calculations are carried out based on the parame¬ 
ters summarized in Tables 1 and 2 which can be derived from 
Refs. [48,93-95,97]. These parameters are kept constant un¬ 
less otherwise mentioned specifically. 

Combining Eqs. (24)-(27), (30), and (33), one can obtain the 
curves of T l v , T 3 )) , and T 4l) varying with the current density i 
for different given values of y, as shown in Fig. 3. When the 
current density is small, Ti_, and T 4> , approach to T 0 and T, 
respectively. When the current density increases, the effi¬ 
ciency of the fuel cell decreases and the waste heat produced 
in the fuel cell increases. In order to efficiently utilize the 
waste heat, (T lii; - T 0 ) and (T - T 4 ,) are required to increase 
with the increase of the current density, and T 3v decreases 
with the decrease of T 4> , accordingto Eq. (21). Substituting T 1)V , 
T 3iV , and T 4> , into Eq. (24), one can obtain the maximum effi¬ 
ciencies ij h m of various heat engines for a given heat input. 
Combining Eqs. (31)—(33), one can generate the curves of the 
power output and efficiency of SOFC-based hybrid systems 
varying with the current density, as shown in Fig. 4, where 
P* = P/A c and P* = P e /A c [96], It can be seen from Fig. 4 that the 
maximum power density of the SOFC is 6622 W nT 2 , which is 


Table 1 — Parameters used in the SOFC model. 

Parameter 

Value 

Molar enthalpy change at 

-249.2 

T = 1223 K, Ah 0 (kj mol -1 ) 


Molar Gibbs free energy 

-180.3 

change at T = 1223 K, Ag° (kj mol -1 ) 


Operating pressure, p 0 (atm) 

1 

Fuel composition pressures, 

0.66; 0.34 

Ph 2 ; Ph 2 o (atm) 


Air composition pressures, 

0.21; 0.79 

Po 2 ;Pn 2 (atm) 


Number of electrons, n e 

2 

Faraday constant, F (C mol -1 ) 

96,485 

Universal gas constant, 

8.314 

R (J mol -1 K -1 ) 


Activation overpotential: 


Anode exchange current 

6610 

density, i 0 , a (A m -2 ) 


Cathode exchange current 

3550 

density, io, c (A m -2 ) 


Charge transfer coefficient 

0.5 

of the electrodes, a 


Ohm overpotential: 


Flowing length in anode: 

0.01; 0.00298; 

5 a (cm); a a (Q cm); b a (K) 

-1392 

Flowing length in cathode: 

0.22; 0.00811; 

5 C (cm); a c (Q cm); b c (K) 

600 

Flowing length in electrolyte: 

0.004; 0.00294; 

<5 e (cm); a e (Q cm); b e (K) 

10,350 

Flowing length in interconnector: 

0.0085; 0.12; 

<5j (cm); a; (Q cm); b; (K) 

4690 

Concentration overpotential: 


Limiting current density of 

2.99 x 10 4 

H 2 , il,h 2 (a irr 2 ) 


Limiting current density of 

2.16 x 10 4 

0 2 , i t ,0 2 (A irr 2 ) 



situated between the experimental values of 
5710-8830 W mr 2 [98,99], Similarly, using Eqs. (24)—(27), 
(30)-(32), and (34), one can generate the curves of the power 
output and efficiency of MCFC-based hybrid systems varying 
with the current density, as shown in Fig. 5. It can be seen 
from Fig. 5 that the maximum efficiency of the MCFC is 0.74, 
which is situated between the experimental values of 0.4-0.86 
[ 100 ], 

It can be seen from Figs. 4 and 5 that the performance of 
high-temperature fuel cells can be greatly improved by 
combining heat engines. For a high-temperature fuel-cell 
heat-engine hybrid system, there always exist one maximum 
power output P max and one maximum efficiency r] max . When 
i > ip, both the power output and the efficiency of the hybrid 
system decrease as the current density is increased. When 
i < i n , both the power output and the efficiency of the hybrid 
system decrease as the current density is decreased. It shows 
that the current density i v at the maximum efficiency and 
current density i P at the maximum power output are, 
respectively, the lower and upper bounds of the optimized 
current density of the high-temperature fuel cell. Thus, the 
optimal region of the current density of the high-temperature 
fuel cell may not be arbitrarily chosen and should be situated 
in the range of i v < i < i P . Consequently, the optimum regions 
of the power output and efficiency of the hybrid system 
should be, respectively, situated in the range of P, < P < P max 
and vp< V < Vmax, where P, and rj P are, respectively, the power 
output at the maximum efficiency and efficiency at the 
maximum power output of the hybrid system. 


4. Several interesting cases 

Discussions above only cover a few examples of hybrid sys¬ 
tems. Several interesting cases of hybrid systems need to be 
further researched. As described above, the unified model 
proposed here also includes the hybrid systems composed of 


Table 2 - Parameters used in the MCFC model. 


Parameter Value 

Molar enthalpy change at 247.4 

T = 923 K, Ah 0 (kj mol" 1 ) 

Molar Gibbs free energy change 197.0 

at T = 923 K, Ag° (kj mol -1 ) 

Partial pressure of H 2 , Ph 2 (atm) 0.60 

Partial pressure of CO 2 in the 0.15 

anode, p C o 2 ,an (atm) 

Partial pressure of CO 2 in the 0.08 

cathode, p C o 2 ,ca (atm) 

Partial pressure of H 2 0, pn 2 o (atm) 0.25 

Partial pressure of 0 2 , po 2 (atm) 0.08 

Activation overpotential: 

Anode exchange current 2404 

density, i 0>a (Am -2 ) 

Cathode exchange current 602 

density, i 0 , c (A m -2 ) 

Ohm overpotential: 

Flowing length in electrolyte: 0.17; 2.94; 

5 e (cm); a e (Q cm); b e (K) 3016 

Reference temperature, T 0 (K) 923 
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(a) 



(b) 



(c) 


Fig. 3 - The curves of (a) T lpn , (b) T 3i „ and (c) T 4j , versus the 
current density of SOFC-based hybrid systems for the 
parameters T = 1223 K, T 0 = 300 K, K„ = K R = K t , I = 0.90, 
x = 0.85, b 3 = 0.50 x 10“ 2 , and b 2 = 0.50 x 10 4 . Curves 1, 
2, and 3 correspond to the cases of (1) y = 1.4, (2) y = 1.0, 
and (3) y = 1/1.4, respectively. 


high-temperature fuel cells and the Carnot cycle, gas turbine, 
Braysson cycle, etc. Furthermore, the unified model can be 
used to discuss the effects of irreversible losses on the per¬ 
formance of hybrid systems. 

4.1. The case of y = 0 

When c L -> <*> while c H remains finite, it is noted that 7 = 0, and 
the cycle consists of a polytropic, an isothermal, and two 


adiabatic processes, constituting exactly the Braysson cycle 
[79], In such a case, Tq = T 6 = T 5 , there does not exist a 
regenerative process, and consequently, T 2 = T 3 . Using the 
similar method, one can obtain the efficiency and heat input 
of the heat engine as 


T 1 ln(T 4 x/T 1 ) 
I'(T 4 - Tq/x) 


and 


(35) 


= _ K H A(T 4 —Tq/x) _ 

qi ln[(T-T 1 /x)/(T-T 4 )] + (K„/K L )T 1 ln(T 4 x/T 1 )/[I'(T 1 -To)]’ 

(36) 

where r = jfic H ln(xT 4 /Tq)/As can be another parameter to 
describe the internal irreversible losses of the working sub¬ 
stance, and As is the rate of the entropy change during the 
low-temperature isothermal process. By searching for the 
optimum values T lpV and T 4> , of Tq and T 4 for a given heat input 
qi, one can obtain the maximum efficiency ijh.m of the heat 



(a) 



(b) 


Fig. 4 - The curves of (a) power density and (b) efficiency 
versus current density of SOFC-based hybrid systems, 
where i P and i,, are the current densities at the maximum 
power density P max and maximum efficiency r/ max , 
respectively. Curves 1-3 correspond to the cases of y = 1.4, 
y = 1.0, and y = 1/1.4, respectively. Curve 4 corresponds to 
the case of the SOFC. The values of other parameters are 
the same as those used in Fig. 3. 
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(b) 

Fig. 5 - The curves of (a) power density and (b) efficiency 
versus current density of MCFC-based hybrid systems, 
where T = 923 K, T 0 = 300 K, K H = K R = K L , I = 0.91, 
x = 0.85, bi = b 2 = 0.01, and i P and i n are the current 
densities at the maximum power density P max and 
maximum efficiency r/ max , respectively. Curves 1-3 
correspond to the cases of y = 1.4, y = 1.0, and y = 1/1.4, 
respectively. Curve 4 corresponds to the case of the MCFC. 


engine. Combining Eqs. (31)—(34), one can obtain the perfor¬ 
mance characteristic curves of high-temperature fuel-cell 
Braysson engine hybrid systems, as shown by curve 5 in Figs. 6 
and 7. 

When I' = 1, the cycle mode of the heat engine is endor- 
eversible and the performance characteristic curves of high- 
temperature fuel-cell endoreversible Braysson engine hybrid 
systems are shown by curve 4 in Figs. 6 and 7. Curves 4 and 5 in 
Figs. 6 and 7 clearly show that the performance of a fuel-cell 
endoreversible Braysson engine hybrid system is much bet¬ 
ter than that of a fuel-cell irreversible Braysson engine hybrid 
system. It is seen from curves 4 and 5 in Figs. 6 and 7 that for a 


r| 

Fig. 6 — The power output versus efficiency curves of SOFC- 
based hybrid systems, where x = 0.50 and I' = 0.85. Curves 
1, 2, 3, 4, 5, 6, and 7 correspond to the hybrid systems 
composed of the high-temperature SOFC and the different 
engines, which are, respectively, (1) Reversible Carnot 
engine, (2) Endoreversible Carnot engine, (3) Irreversible 
Carnot engine, (4) Endoreversible Braysson engine, (5) 
Irreversible Braysson engine, (6) Endoreversible Ericston 
engine, and (7) Irreversible Ericston engine. The values of 
other parameters are the same as those used in Fig. 3. 


engine hybrid system increases, and the hybrid system may 
be operated around the maximum power output to achieve 
the optimal performance; while the MCFC Braysson engine 
hybrid system may be operated in the range from the power 
output at the maximum efficiency to the maximum power 
output to achieve the optimal performance. 

4.2. The case of y —> °° 


When c H -> oo while c L remains finite, we note that y -> , 

T 2 = T 3 = T 4 , and consequently, T 6 = T 5 . The heat engine cycle 
consists of an isothermal, an isobaric, and two adiabatic 
processes. It is based on an Ericsson cycle for the high- 
temperature heat addition and a Brayton cycle for the low- 
temperature heat rejection. Although such a cycle has not 
been named, its performance has been evaluated [81]. Similar 
to the Braysson cycle, such a cycle may be simply referred to 
as the Ericston cycle. Using the similar method, one can obtain 
the efficiency and heat input of the heat engine as 


x(T 5 -TQ 

,h T 1 I'ln(T 5 /T 1 ) 

and 


(37) 


_KhA_ 

(1/(T - Ti/x)) + (K H /K L )xln[(T 5 - T 0 )/(T t - ToMI'TjnfTs/Ta)] ’ 


(38) 


SOFC-based hybrid system, the power output will increase where I' = As/[mc L ln(T 5 /Ti)], As is the rate of the entropy 
and then decrease as the efficiency of the SOFC Braysson change during the high-temperature isothermal process. 
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Fig. 7 - The power output versus efficiency curves of 
MCFC-based hybrid systems where x = 0.50 and I' = 0.85. 
Curves 1, 2, 3, 4, 5, 6, and 7 correspond to the hybrid 
systems composed of the high-temperature MCFC and the 
different engines, which are, respectively, (1) Reversible 
Carnot engine, (2) Endoreversible Carnot engine, (3) 
Irreversible Carnot engine, (4) Endoreversible Braysson 
engine, (5) Irreversible Braysson engine, (6) Endoreversible 
Ericston engine, and (7) Irreversible Ericston engine. The 
values of other parameters are the same as those used in 
Fig. 5. 


Combining Eqs. (31)—(34), (37), and (38), one can obtain the 
performance characteristic curves of high-temperature fuel¬ 
cell Ericston heat engine hybrid systems, as shown by curve 7 
in Figs. 6 and 7. Similarly, the performance characteristic 
curves of high-temperature fuel-cell endoreversible Ericston 
engine hybrid systems can be plotted by using the condition 
I' = 1, as shown by curve 6 in Figs. 6 and 7. Curves 6 and 7 in 
Figs. 6 and 7 clearly show that the fuel-cell Ericston engine 
hybrid system may be operated around the maximum power 
output to achieve the optimal performance. 


4.3. The case of c H -*• °° and c L -» °° 


When c H —► 00 and c L —> °=, it follows that Ti = T 6 = T 5 , 
T 2 = T 3 = T 4 , and the heat engine is a Carnot engine. One can 
obtain the efficiency and heat input of the Carnot engine as [85] 



and 


(39) 


_ K h A(T - T„) _ 

1 + {(Kh/K l )T!(T - T 4 )/[I'T 4 (T 1 - T 0 )]} ’ 


(40) 


where I' is the ratio of entropy change rates during high- and 
low-temperature isothermal processes. Combining Eqs. 
(31)-(34), (39), and (40), one can obtain the performance 
characteristic curves of high-temperature fuel-cell Carnot 
engine hybrid systems, as shown by curve 3 in Figs. 6 and 7. By 
substituting I' = 1 into Eqs. (39) and (40), we can also evaluate 
the performance of high-temperature fuel-cell endoreversible 
Carnot engine [86,101] hybrid systems and performance 
characteristic curves are shown by curve 2 in Figs. 6 and 7. 


These characteristics are also suitable for the hybrid systems 
composed of the high temperature fuel cell and the endor¬ 
eversible Stirling or Ericsson cycle with the ideal regeneration 
[72], If we further assume that two heat transfer processes 
between the working substance in the Carnot engine and two 
reservoirs are ideal, i.e., K H = K L —* °°, Ti = T 0 , and T = T 4 , the 
cycle of the heat engine is reversible [28,102], the Carnot effi¬ 
ciency is entirely determined by the temperatures of the fuel 
cell and environment as 

4h = l-y- (41) 

Combining Eqs. (31)—(34), and (41), one can obtain the 
performance characteristic curves of high-temperature fuel¬ 
cell reversible Carnot engine hybrid systems, as shown by 
curve 1 in Figs. 6 and 7. These characteristics are also suitable 
for the hybrid systems composed of the high temperature fuel 
cell and the reversible Stirling or Ericsson cycle with the ideal 
regeneration [72]. It can be seen from curves 1, 2, and 3 in Figs. 
6 and 7 that the performance of a fuel-cell ideal Carnot engine 
hybrid system is much better than that of a fuel-cell irre¬ 
versible Carnot engine hybrid system. 


4.4. The case of y = 1 


When c H = c L = c p or c H = c L = c„, it is noted that y = 1 and the 
cycle mode of the heat engine may be the Brayton or Otto 
cycle. For a Brayton cycle, the rate of heat flow in three pro¬ 
cesses can be expressed as [80] 


K H A H (T 4 - T 3 ) 
ln[(T - T 3 )/(T - T 4 )] 


m Cp(T 4 -T 3 ), 


(42) 


_ K l A l (T 6 - TQ 
H2 ln[(T 6 - T 0 )/(Ti - T 0 )] 


bic p (T 6 - T a ), 


and 


(43) 


<5r — KrA.r(T 5 — T 3 ) — K R A R (T 6 — T 2 ) = mc p (T 3 — T 2 ) 

= mc p ( T 5 -T 6 ), (44) 

respectively, and the general performance characteristics of 
such hybrid systems have been shown by curve 2 in Figs. 4 and 
5. 

4.4.1. The case of 1 = 1 

When I = 1, the Brayton cycle with the regenerative process is 
endoreversible, K R —> °°, T 5 = T 3 , T 6 = T 2 , and T 2 /T 4 = TVT 5 . The 
efficiency and heat input rate of the endoreversible Brayton 
cycle are, respectively, 

^ = ( 45 > 

and 

= _ AK h T 4 (1 - x) _ 

qi ln[(T - xT 4 )/(T - T 4 )] + (K H /K L )ln[(T a /x - To)/(T 4 - T 0 )]' 

(46) 

Combining Eqs. (31)—(34), (45), and (46), one can obtain the 
performance characteristic curves of high-temperature fuel- 
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Fig. 8 - The power output versus efficiency curves of SOFC- 
Brayton engine hybrid systems under different conditions, 
where x = 0.80. Curves 1, 2, 3, 4, 5, 6, and 7 correspond to 
the cases of (1) Endoreversible Brayton cycle, (2) 
Endoreversible gas turbine, (3) Irreversible gas turbine, (4) 
Endoreversible Brayton cycle without the regenerative 
process, (5) Irreversible Brayton cycle without the 
regenerative process, (6) Endoreversible gas turbine 
without the regenerative process, and (7) Irreversible gas 
turbine without the regenerative process. The values of 
other parameters are the same as those used in Fig. 3. 


cell endoreversible Brayton engine hybrid systems, as shown 
by curve 1 in Figs. 8 and 9. 

4.4.2. The case without regenerative process 
If there does not exist a regenerator in the Brayton cycle jl03], 
K r = 0, T 2 = T 3 , T 6 = T 5 , and the waste gas from the turbine is 
directly released to the environment. In such a case, I = (TV 
T 5 )(T 4 /T 2 ) and one can obtain the efficiency and heat input rate 
of the cycle as 



Fig. 9 - The power output versus efficiency curves of 
MCFC-Brayton engine hybrid systems under different 
conditions, where x = 0.80. Curves 1, 2, 3, 4, 5, 6, and 7 
correspond to the cases of (1) Endoreversible Brayton cycle, 
(2) Endoreversible gas turbine, (3) Irreversible gas turbine, 
(4) Endoreversible Brayton cycle without the regenerative 
process, (5) Irreversible Brayton cycle without the 
regenerative process, (6) Endoreversible gas turbine 
without the regenerative process, and (7) Irreversible gas 
turbine without the regenerative process. The values of 
other parameters are the same as those used in Fig. 5. 


endoreversible Brayton engine [104] hybrid systems and 
obtain its performance characteristic curves, as shown by 
curve 4 in Figs. 8 and 9. 

Comparing curves 1 and 4 in Figs. 8 and 9, one can find that 
the introduction of a regenerator can effectively utilize the 
waste heat from the turbine and improve the performance of 
the hybrid system [80]. 

It should be pointed out that, as long as c p is replaced by c„, 
the results obtained in Sections 4.4—4.4.2 are also suitable for 
the hybrid systems composed of the high temperature fuel 
cell and the Otto cycle. 


Vh 


x((T 4 /I) - (Vx)) 
T 4 - Ti/x 


and 


(47) 


AK H (T 4 — T 4 /x) 

qi ln[(T — Ti/x)/(T — T 4 )[ + (K„/K L )ln[(xT 4 /I - T 0 )/(T! - T 0 )[' 

(48) 

Combining Eqs. (31)—(34), (47), and (48), one can obtain the 
performance characteristic curves of high-temperature fuel- 


4.4.3. The case of T 2 = T 0 

When Ti = T 0 and the cycle is an open Brayton cycle, the cycle 
model can be used to discuss the performance of high- 
temperature fuel-cell gas-turbine hybrid systems 105,106]. 
The efficiency and heat input rate of the gas turbine cycle are 


Vh = 


T 0 (T 4 /I - T 3 ) 
T 3 (T 4 -T 3 ) 


and 


(49) 


= _ AK H (T 4 - T 3 ) _ 

qi ln[(T - T 3 )/(T - T 4 )] + (K h /Kr)(T 3 - T 0 /x)/[T 0 T 4 /(IT 3 ) - (T 0 /x)] ’ 


(50) 


cell Brayton engine hybrid systems, as shown by curves 4 and 
5 in Figs. 8 and 9. Substituting I = 1 into Eqs. (47) and (48), one 
can discuss the performance of the high-temperature fuel-cell 


respectively, and the performance characteristic curves of 
hybrid systems are shown by curve 3 in Figs. 8 and 9. Similarly, 
one can also obtain the performance characteristic curves of 
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high-temperature fuel-cell endoreversible gas turbine hybrid 
systems by introducing the conditions K R -» °°, T 5 = T 3 , T 6 = T 2 , 
and T 2 /T 0 = T 4 /T 5 , as shown by curve 2 in Figs. 8 and 9. 
Comparing curves 1,2, and 3 in Figs. 8 and 9, one can obtain that 
the performance of a fuel-cell endoreversible gas turbine hybrid 
system is not only much better than that of a fuel-cell irrevers¬ 
ible gas turbine hybrid system but also slightly better than that 
of a fuel-cell endoreversible Brayton engine hybrid system. 

If there does not exist a regenerator in the gas turbine 107], 
one can obtain the efficiency and heat input rate of the cycle 
as 

.. , x((T 4 /I) - (T 0 /x)) 

Vh = 1 -T 4 — To/x (51) 


and 

AK h (T 4 - Tq/x) 
qi ln[(T - T 0 /x)/(T - T 4 )] ’ 


(52) 


respectively. The performance characteristic curves of the 
hybrid systems composed of the high-temperature fuel cell 
and the gas turbine without the regenerative process are 
shown by curve 7 in Figs. 8 and 9, and the performance 
characteristic curves of high-temperature fuel-cell endor¬ 
eversible gas turbine hybrid systems are shown by curve 6 in 
Figs. 8 and 9. Curves 4-7 in Figs. 8 and 9 show that if there is 
not any regenerator in the engine cycle, the performance of 
the fuel-cell gas turbine hybrid system is close to that of a 
fuel-cell Brayton engine hybrid system. 



Fig. 10 — The power output versus efficiency curves of 
high-temperature fuel-cell Diesel engine hybrid systems 
under different conditions, where x = 0.70. Curves 1, 2, and 
3 are the performance curves of SOFC-Diesel engine hybrid 
systems which correspond to the cases of (1) 
Endoreversible Diesel cycle, (2) Endoreversible Diesel cycle 
without the regenerative process, and (3) Irreversible 
Diesel cycle without the regenerative process. Curves T, 2', 
and 3' are the performance curves of the MCFC-Diesel 
engine hybrid system which correspond to the same cases 
as curves 1, 2, and 3, respectively. The values of other 
parameters are the same as those used in Figs. 3 and 5. 


4.5. The case of y = Cp/c v 


When c H = c p and c L = c„, clearly, y = 1.4 for diatomic gases 
[108]. In such a case, the hybrid system is composed of the 
high-temperature fuel cell and the Diesel engine. The general 
performance characteristic curves of the hybrid system have 
been shown by curve 1 in Figs. 4 and 5. Substituting y = 1.4, 
I = 1, and (T 3 /T 2) 1 ' 4 = Ts/T 6 into Eqs. (24) and (25), one can 
obtain the efficiency and heat input rate of the endoreversible 
Diesel engine as 


Ti[(T;/T3) 14 -l] 
1-4(T 4 - T 3 ) 


(53) 


and 

= _ AK h (T'-T 3 ) _ 

^ i- T ~T 3 . K » ^(xn/TQ 14 -^ , K h ^(n/T,) 14 —To ’ 

T—T 4 0.4K r Tl (r 4 /T 3 ) 1 - 4 -T 1 /x + 1-4Kl T 4 -T 0 

(54) 


where r 4 =[(1.4((T 3 /T 1 )-l/x))/((T 1 /x)- 1 - 4 -T^ 1 - 4 ))] 1/1 - 4 . Com¬ 
bining Eqs. (31)—(34), (53), and (54), one can obtain the per¬ 
formance characteristic curves of high-temperature fuel-cell 
endoreversible Diesel engine hybrid systems, as shown by 
curves 1 and T in Fig. 10. Curves 1 and T in Fig. 10 and curve 1 
in Figs. 4(a) and 5(a) show that the decrease of x may cause the 
reduction of the power output of a fuel-cell Diesel engine 
hybrid system. 



Fig. 11 - The power output versus efficiency curves of 
high-temperature fuel-cell Atkinson engine hybrid 
systems under different conditions, where x = 0.70. 
Curves 1, 2, and 3 are the performance curves of the SOFC- 
Atkinson engine hybrid system which correspond to the 
cases of (1) Endoreversible Atkinson cycle, (2) 
Endoreversible Atkinson cycle without the regenerative 
process, and (3) Irreversible Atkinson cycle without the 
regenerative process. Curves 1', 2', and 3' are the 
performance curves of the MCFC-Atkinson engine hybrid 
system which correspond to the same cases as curves 1, 2, 
and 3, respectively. The values of other parameters are the 
same as those used in Figs. 3 and 5. 
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If there does not exist a regenerator in the Diesel cycle, 
K r = 0, T 2 = T 3 , T 6 = T s , and I = (T 1 /T 5 )(T 4 /T 2 ) 1 - 4 . From Eqs. 
(21)-(23), one can obtain the efficiency and heat input rate of 
the cycle as 


If there does not exist a regenerator in the Atkinson cycle 
[82], then K R = 0, T 2 = T 3 , T 6 = T 5 , and I = (T 1 /T 5 )(T 4 /T 2 ) V7 . From 
Eqs. (21)— (23), one can obtain the efficiency and heat input 
rate of the cycle as 


T 1 (xT 4 /T 1 ) 1 ' 4 /I-T 1 
Vh 1.4(T 4 - Tj/x) 

and 


( 55 ) 


Vh = 1 


7Ti ^(xT 4 /T 1 ) 5/7 /I — 1 
5(T 4 - Tj/x) 


(59) 


_AK h (T 4 - Tt/x)_ 

ln[(T - T!/x)/(T - T 4 )] + [K„/(1.4K L )]ln{ [(T 1 (xT 4 /T 1 ) 1 - 4 )/l - T 0 ] /(T 4 - T 0 )} ' 


(56) 


and 


_ AK h (T 4 - T t /x) _ 

ln[(T - T 4 /x)/(T - T 4 )] + [7K H /(5K L )]ln{ [(t^xTi/Tj) 577 )/! - T 0 ] /(T 4 - T 0 )} ' 


(60) 


Using Eqs. (31)—(34), (55) and (56), one can obtain the per¬ 
formance characteristic curves of high-temperature fuel-cell 
Diesel engine hybrid systems, as shown by curves 3 and 3' in 
Fig. 10. Substituting I = 1 into Eqs. (55) and (56), one can obtain 
the hybrid systems composed of the high-temperature fuel 
cell and the endoreversible Diesel engine with no regenerator 
[109], and performance characteristic curves are shown by 
curves 2 and 2' in Fig. 10. Curves 1, T, 2, and 2' in Fig. 10 show 
the excellence of the fuel-cell Diesel engine hybrid systems 
with the regenerator. 


4.6. The case of y = c v /c p 


When c H = c„ and c L = c p , we note that y = 5/7 for diatomic 
gases. In such a case, the hybrid system is composed of the 
high-temperature fuel cell and the Atkinson engine. The 
general performance characteristic curves of hybrid systems 
have been shown by curve 3 in Figs. 4 and 5. Substituting y = 5/ 
7, I = 1, and (T 3 /T 2 ) s/7 = T 5 /T 6 into Eqs. (24) and (25), one can 
obtain the efficiency and heat input rate of the endoreversible 
Atkinson engine as 


Vh 


7t 1 [(tvt 3 ) 5/7 -i] 

50F 4 -T 3 ) 


(57) 


and 


= _ AK h (T 4 —T 3 ) _ 

^ , T —T 3 7K„ T 1 (T' i /T 3 ) s/7 -T 1 /x 7K H T 4 (T 4 /T 3 ) 5/7 —T 0 ’ 

T-T' + 2K r Tl (xTi/T 4 ) 5/7 -T 3 5K l T 4 -T 0 

(58) 

where T t = [(((Ts/T,) - (l/x))(5/7))/((T a /x)- s/7 -T 3 5/7 )] 7/5 . 

Combining Eqs. (31)—(34), (57), and (58), one can obtain the 
performance characteristic curves of high-temperature fuel¬ 
cell endoreversible Atkinson engine hybrid systems, as shown 
by curves 1 and T in Fig. 11. 


Using Eqs. (31)—(34), (59), and (60), one can obtain the per¬ 
formance characteristic curves of high-temperature fuel-cell 
Atkinson engine hybrid systems, as shown by curves 3 and 3' 
in Fig. 11. Substituting I = 1 into Eqs. (59) and (60), one can 
discuss the performance of the hybrid systems composed of 
the high-temperature fuel cell and the endoreversible Atkin¬ 
son engine without a regenerative process and can generate 
the performance characteristic curves, as shown by curves 2 
and 2' in Fig. 11. The curves in Fig. 11 show that the power 
output will increase and then decrease as the efficiency of a 
fuel-cell Atkinson engine hybrid system increases, and the 
hybrid system may be operated around the maximum power 
output to achieve the optimal performance of the hybrid 
system. 


5. Conclusions 

In this study, the performances of various high-temperature 
fuel cells, heat engines, and their hybrid systems are 
reviewed. A new unified model of high-temperature fuel-cell 
heat-enigne hybrid systems is proposed. In detail, the per¬ 
formance of hybrid systems is studied by considering the ir¬ 
reversibilities primarily caused by overpotentials and heat 
leakage in high-temperature fuel cells, finite-rate heat trans¬ 
fer between the working substance and the heat reservoirs, 
and internal losses of the working substance in heat engines. 
The performance characteristic curves of high-temperature 
fuel cells and hybrid systems are, respectively, obtained. The 
optimally working regions of some of the important perfor¬ 
mance parameters of hybrid systems are determined. Some 
interesting cases of hybrid systems are also discussed in 
detail. The results obtained show that hybrid systems can 
effectively utilize the high temperature waste heat and 
improve both the power output and the efficiency of high- 
temperature fuel cells. 
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